Lipidomic studies revealing serological markers associated with the occurrence of retinopathy in type 2 diabetes

Purpose The duration of type 2 diabetes mellitus (T2DM) and blood glucose levels have a significant impact on the development of T2DM complications. However, currently known risk factors are not good predictors of the onset or progression of diabetic retinopathy (DR). Therefore, we aimed to investigate the differences in the serum lipid composition in patients with T2DM, without and with DR, and search for potential serological indicators associated with the development of DR. Methods A total of 622 patients with T2DM hospitalized in the Department of Endocrinology of the First Affiliated Hospital of Xi’an JiaoTong University were selected as the discovery set. One-to-one case–control matching was performed according to the traditional risk factors for DR (i.e., age, duration of diabetes, HbA1c level, and hypertension). All cases with comorbid chronic kidney disease were excluded to eliminate confounding factors. A total of 42 pairs were successfully matched. T2DM patients with DR (DR group) were the case group, and T2DM patients without DR (NDR group) served as control subjects. Ultra-performance liquid chromatography–mass spectrometry (LC–MS/MS) was used for untargeted lipidomics analysis on serum, and a partial least squares discriminant analysis (PLS-DA) model was established to screen differential lipid molecules based on variable importance in the projection (VIP) > 1. An additional 531 T2DM patients were selected as the validation set. Next, 1:1 propensity score matching (PSM) was performed for the traditional risk factors for DR, and a combined 95 pairings in the NDR and DR groups were successfully matched. The screened differential lipid molecules were validated by multiple reaction monitoring (MRM) quantification based on mass spectrometry. Results The discovery set showed no differences in traditional risk factors associated with the development of DR (i.e., age, disease duration, HbA1c, blood pressure, and glomerular filtration rate). In the DR group compared with the NDR group, the levels of three ceramides (Cer) and seven sphingomyelins (SM) were significantly lower, and one phosphatidylcholine (PC), two lysophosphatidylcholines (LPC), and two SMs were significantly higher. Furthermore, evaluation of these 15 differential lipid molecules in the validation sample set showed that three Cer and SM(d18:1/24:1) molecules were substantially lower in the DR group. After excluding other confounding factors (e.g., sex, BMI, lipid-lowering drug therapy, and lipid levels), multifactorial logistic regression analysis revealed that a lower abundance of two ceramides, i.e., Cer(d18:0/22:0) and Cer(d18:0/24:0), was an independent risk factor for the occurrence of DR in T2DM patients. Conclusion Disturbances in lipid metabolism are closely associated with the occurrence of DR in patients with T2DM, especially in ceramides. Our study revealed for the first time that Cer(d18:0/22:0) and Cer(d18:0/24:0) might be potential serological markers for the diagnosis of DR occurrence in T2DM patients, providing new ideas for the early diagnosis of DR. Supplementary Information The online version contains supplementary material available at 10.1186/s12967-024-05274-9.


Introduction
Type 2 diabetes mellitus (T2DM) is a common chronic disease in many countries, and its prevalence is growing as people's lifestyles are changing [1].Diabetes causes various complications, classified as either macrovascular complications (such as cardiovascular disease and stroke) or microvascular complications (such as kidney disease) [2].Diabetic retinopathy (DR), a specific microvascular complication of diabetes, is the most common cause of vision loss in people of working age [3,4].Poor glycemic control, hypertension, and diabetes duration are major risk factors for DR [5].Although intensive risk factor control reduces the risk of DR progression and vision loss, many diabetic patients continue to develop DR with strict glycemic and blood pressure control [6].Despite increasing research supporting the efficacy of routine DR screening to prevent DR and early treatment to reduce the risk of vision loss, there are no specific biomarkers for diagnosing the onset and early progression of DR.Additionally, new and more effective strategies are awaited to prevent and treat the progression of DR.
Accumulating evidence suggests that disruption in lipid metabolism is an early event in the pathogenesis of diabetes complications.Previous studies found that levels of multiple lipid species, including glycerophospholipids, sphingolipids and glycerolipids, are critical risk factors for T2DM and its complications [7,8].Lysophosphatidylcholine (LPC) is a main glycerophospholipid known for its essential role in lipid and glucose metabolism, and LPC has been intensively studied in the development of metabolic diseases including T2DM [9].Sphingolipids, including ceramides (Cer), sphingomyelins (SM) and gangliosides, have a variety of intra-and extracellular effects on glucose homeostasis and metabolic disease [10] Numerous studies suggest Cer, a crucial lipid intermediate in sphingolipid metabolism, is a major contributing factor for insulin resistance, and inhibition or depletion of enzymes driving de novo ceramide synthesis can prevent the development of diabetes in mice [7,11,12].In contrast, a decrease in very long chain Cer is correlated with the development of macroalbuminuria in diabetes [13].Accelerated sphingolipid catabolism' leading to an increase in glucosylceramide or glycosphingolipids might contribute to the neuronal pathologies of DR [14].In addition, SM produced by the transfer of a phosphocholine moiety from phosphatidylcholine to the ceramide backbone has been linked to insulin resistance [15,16] and is also an independent marker of cardiovascular disease [17].Thus, dysregulated lipid metabolism is a major contributor to the pathogenesis of T2DM and its complications, and specific lipid species that are responsible for the occurrence of DR are rather obscure.
Lipidomics offers solid platforms for identifying novel lipid mediates in biochemical processes of lipid metabolism, thus providing new opportunities for disease prediction and detection [18,19].Lipidome analysis is performed by liquid chromatography and electrospray ionization-tandem mass spectrometry (LC-MS/MS) for molecular lipid identification and quantification and multiple reaction monitoring (MRM) for targeted quantification of those lipid species.Lipid-based biomarkers offer unique options for precision medicine by providing sensitive diagnostic tools for disease prediction and monitoring [20].Using a quantitative metabolomics approach, Emil et al. compared the aqueous humor and serum concentrations of metabolites in senior adults with an without diabetes who underwent cataract surgery [21].However, the field of lipidomics studies of DR is still in its early stages, with few studies published and little replication of results [22].
In this study, we aimed to find reliable serum lipidbased biomarkers for the presence of DR in patients with T2DM by using two cohorts.To this end, serum samples of the discovery cohort was subjected to untargeted lipidomics analysis to search for differentially abundant lipids between individuals without and with DR.In the validation cohort, the observed differential lipid molecules were validated using mass spectrometry MRM targeting techniques.We hypothesized that DR has a distinctive serum lipid signature and that particular lipid species can act as biomarkers for T2DM patients with DR.

Participants
A total of 622 participants with T2DM hospitalized in the Endocrinology Department of the First Affiliated Hospital of Xi'an JiaoTong University were screened as the discovery set.Participants with chronic kidney disease [estimated glomerular filtration rate (eGFR) < 90 (mL/min/1.73m 2 )] were excluded from the selection.We conducted pair matching according to the traditional risk factors for DR (including age, duration of diabetes, HbA1c level, and hypertension).For the discovery cohort, we selected 42 T2DM patients with DR (DR group).The control participants were 42 T2DM patients without DR (NDR group), and they were matched to patients in the DR group by age (in 5-year bands), diabetes duration (in 5-year bands), HbA1c levels (in 0.5% bands), and hypertension status.
Lipid markers of DR identified from the discovery cohort were quantified in a separate sample cohort (validation cohort).We first screened 531 T2DM patients.Individuals with chronic kidney disease [eGFR < 90 (mL/ min/1.73m 2 )] were excluded from the selection.Then, we conducted 1:1 propensity score matching (PSM) (matching tolerance = 0.02) by age, diabetes duration, HbA1c level, hypertension status, sex, BMI, systolic blood pressure (SBP), diastolic blood pressure (DBP), and eGFR.For the validation cohort, 95 T2DM patients with DR (DR group) and 95 T2DM patients without DR (NDR group) were included.

Lipid extraction
The serum samples were thawed slowly at 4 °C, 100 µL of the sample was placed in a 96-well plate, 300 µL of isopropanol (prechilled at -20 °C) spiked with internal standards (SPLASH® LIPIDOMIX® Mass Spec Standard, Avanti, USA) was added, and the samples were vortexed and mixed for 1 min and then centrifuged at 4 °C for 20 min at 4000 rcf after resting overnight at -20 °C as previously reported [23].The supernatant was injected for LC-MS/MS analysis, and 10 µL of each supernatant was mixed into quality control (QC) samples to assess the reproducibility and stability of the LC-MS analysis process.
All samples were analyzed in data-dependent acquisition (DDA) mode with the following positive/negative ionization settings: spray voltage, 3.8/-3.2kV; aux gas heater temperature, 350 °C; and capillary temperature, 320 °C.The full scan mass range was 200-2000 m/z with 70,000 mass resolution at m/z 200 and AGC set to 3e6 with a maximum ion injection time of 100 ms.The top three precursors were selected for subsequent MS fragmentation with a maximum ion injection time of 50 ms and resolution of 17,500 at m/z 200, and the AGC was 1e5.The stepped normalized collision energy was set to 15, 30, and 45 eV.

Data preprocessing and quality control
The raw data obtained from the LC-MS/MS detection were imported into LipidSearch v.4.1 (Thermo Fisher Scientific, USA) for lipid identification and quantification.The following parameters were used for lipid identification and peak extraction: the type of identification was Product, the mass deviation of the parent and daughter ions was 5 ppm, and the response threshold was set to 5.0% of the relative response deviation of the daughter ions; the quantitative parameters were set to calculate the peak areas of all identified lipids, and the peak extraction mass deviation was set to 5 ppm.For ESI + data, [M + H]+, [M + NH4]+, and [M + Na] + were selected as adducts, while for ESI-data, [M-H]-, [M-2 H]-, and [M-HCOO]were selected as adducts.The peak alignment was performed for all identified lipids, and those not marked as "rejected" were considered for inclusion in the subsequent analysis.
For data preprocessing, raw data exported from Lipid-Search were further analyzed by meta X [25].The data preprocessing included (1) Removing lipid molecules with more than 50% missing information in QC samples and more than 80% missing information in experimental samples (i.e., LipidIon in the table); (2) Filling the missing values using the k-nearest neighbor (KNN) algorithm; (3) Correcting the batch effect using quality control-based robust LOESS signal correction (QC-RLSC); (4) Using probabilistic quotient normalization (PQN) to normalize the data to obtain the relative peak areas; and (5) Removing the lipid molecules with a coefficient of variation (CV) greater than 30% of the relative peak areas from all QC samples.Data quality was assessed by the reproducibility of QC sample assays.The assessment included chromatogram overlap of QC samples, principal component analysis (PCA), number of extracted peaks, and differences in peak response intensity.

Data processing
A combination of multivariate statistical analysis and univariate analysis was used to screen for lipids of which the abundance differed between groups.The multivariate statistical analysis methods used were principal component analysis (PCA) and partial least squares methoddiscriminant analysis (PLS-DA).PCA is an unsupervised pattern recognition method, and PLS-DA is a supervised pattern recognition method.The univariate analyses were fold change (FC) and Student's t test.The FC was obtained by fold change analysis, and the p value pairs of the t test were corrected for the false discovery rate (FDR) to obtain a q-value.The differential lipid molecule screening conditions were as follows: (1) variable importance in the projection (VIP) ≥ 1 for the first two principal components of the PLS-DA model; (2) fold change ≥ 1.2 or ≤ 0.83; and (3) p value < 0.05.

Targeted lipid quantification by MRM in validation samples
The identified differential lipids were further quantified by multiple reaction monitoring (MRM).For lipid extraction, the procedure was consistent with the untargeted experiment as described.The MRM transition list is shown in Table S1.For MRM quantification, all validation samples were analyzed on a QTRAP 5500 mass spectrometer with a CSH C18 column (1.7 μm 2.1*100 mm, Waters, USA) for separation.All lipids were subjected to targeted quantification in ESI + mode with a specific transition setting.

Statistical analysis
The clinical data of samples are presented as the mean ± standard deviation (SD) for normally distributed variables or the median (interquartile range) for abnormal distribution.Comparisons between the case group and the control group were made using a two-tailed t test or Mann-Whitney U test for continuous data and the X 2 test for categorical data.The calculation of the area under the curve (AUC) in receiver operating characteristic (ROC) curve analysis was used to evaluate the discriminatory ability of the markers.Logistic regression models were applied to assess the relationship between lipid molecules and the presence of DR.The odds ratios (ORs) with 95% confidence intervals (CIs) were calculated for the molecules with 1-SD changes.The known risk factors for DR, such as CHOL, TGs, LDL-c, and HDL-c, were added to multivariate logistic regression to calculate the adjusted odds ratios.Ordinal logistic regression models were used to assess the relationships between lipid molecules and DR stages [NDR, nonproliferative DR (NPDR) and proliferative DR (PDR)].

Characteristics of the discovery cohort
Table 1 shows the clinical characteristics of individuals selected for the discovery cohort.There were no significant differences in age and sex between the DR and NDR groups.In fact, these groups were comparable for most metabolic characteristics, such as BMI, diabetes duration, and HbA1c, and there were no significant between-group differences for hypertension status, antihypertensive agent use, hypoglycemic therapy status or NSAID use.The blood pressure and glucose of the participants were treated and controlled.Compared with control subjects with T2DM, T2DM patients with DR had higher levels of LDL-c levels, AST, TBIL, and BUN (Table 1 and Table S2).

Untargeted lipidome-derived biomarkers for diabetic retinopathy: results from the discovery cohort
A total of 1721 lipids were detected.The number of lipids with an RSD (CV) less than or equal to 30% in the QC samples was 1421.The ratio of the number of lipids with CV less than or equal to 30% to the number of all detected lipids in QC samples was 81%.
Fifteen candidate lipids were identified from the discovery cohort.Compared with those of the NDR group, the levels of three Cer and seven SM were significantly lower in the DR group.In contrast, two SM, two LPC and one PC were significantly higher in the DR group (Fig. 1A  and B).More specifically, compared with T2DM patients without DR, T2DM patients with DR showed lower levels of Cer(d18:0/24:0), Cer(d18:0/22:0), Cer(d42:3), SM(d22:0/16:0), SM(d18:1/24:1), SM(d42:0), SM(d40:0), SM(d39:0), SM(d38:0), and SM(d36:0), and higher levels of SM(d20:1/16:1), SM(d34:1), LPC(18:2), LPC(16:0) and PC(34:2).The heat map shows the distribution of these lipids between individuals of the NDR and DR groups (Fig. 1C).The results of ROC analysis and the odds ratios of the lipid markers in the basic logistic regression models are shown in Table 2.The AUC values for the 15 lipids ranged from 0.72 to 0.94.All lipids retained significant ORs after adjusted for CHOL, TG, LDL-c, and HDL-c (adjusted ORs are shown in Table 2).Furthermore, we used ordinal logistic regression, which estimated the odds of being in one higher category of the DR stage (from NDR to PDR) for lipid species, to test the associations between lipid species and DR stage (Table S3; n = 42 in the NDR group, n = 37 in the NPDR group, n = 5 in the PDR group), and we analyzed the data while excluding participants with diabetic macular edema (DME) (n = 4 in the DR group), as before, all lipids retained significant ORs (Table S4).

Characteristics of the validation cohort and targeted lipidomics analysis
The 15 differential lipids found from the discovery cohort were validated in another set of samples.The clinical characteristics of individuals selected for the validation cohort are shown in Table 3.Most metabolic and clinical features were comparable (Table S5), and there was no significant difference in LDL-c between the DR and NDR groups.

Discussion
DR is the most common microvascular complication of diabetes and the main factor contributing to visual impairment in working-age individuals [3].T2DM patients often develop DR despite of proper control of systemic risk factors, indicating the involvement of other pathogenic factors for DR development.To find new and more effective strategies for preventing and treating DR, it is necessary for us to identify novel biomarkers for DR screening or detection.Lipidomics will aid in understanding the mechanism of DR at various stages of the disease, early diagnosis, and the identification of new therapeutic targets.In this study, by using two clinical All data are presented as the mean ± standard deviation (SD) for normally distributed data or the median (interquartile range) for abnormal distribution.Comparisons between groups were made using a two-tailed t test or Mann-Whitney U test for continuous data and the X 2 test for categorical data cohorts, we found that the serum lipidomic profiles in T2DM patients with DR showed significant differences from those in T2DM patients without DR.The differential lipid species in the DR group were linked to disturbances in sphingolipid metabolism.Compared with those in the NDR group, the levels of Cer(d18:0/24:0) and Cer(d18:0/22:0) were significantly lower in the DR group after adjusting for covariates, i.e. known risk factors in both the discovery and validation cohorts.These findings suggest that these two lipid species may be potential serological markers for the diagnosis of DR in patients with T2DM.
In this study, we found two ceramide molecules that were significantly lower in T2DM patients with DR, indicating that they may have disturbed ceramide metabolism compared to T2DM patients without DR.Ceramide is sphingolipid [11] and can be found in VLDL, LDL, and HDL.Consistent with our findings, Fort et al. found a significantly lower abundance of Cer in central retinal tissue obtained postmortem from T2DM patients with DR compared to those without DR [26].Similarly, ceramide levels were shown to be lower and glucosylceramide levels higher in the retinas of diabetic rodents [27].This indicates that diabetes reduces the retinal ceramide content and may suggest that dysregulated sphingolipid metabolism may cause retinal resistance to insulin action [27].These findings imply that ceramide is diverted from the overall pools of retinal sphingolipids toward the glycosylated forms due to hyperglycemia.In contrast, Levitsky et al. found that diabetes-induced increases in mitochondrial ceramide led to impaired mitochondrial function in the retinal pigment epithelial (RPE) cells of the retina [28], and disruption of the blood-retinal barrier might be caused by diabetes-induced overexpression of acid sphingomyelinase.Additionally, inflammation is a common underlying factor in DR, and inflammation generates Cer from SM in the serum membrane.This induces death receptor ligand formation and leads to apoptosis of RPE and photoreceptor cells [29].In addition to diabetes, circulating Cer was shown to strongly correlate with future adverse cardiovascular events.It has recently been discovered that in individuals with atherosclerotic CVD, serum levels of specific Cer species can predict the future risk of cardiovascular death.In the Corogene study, higher concentrations of Cer(d18:1/16:0), Cer(d18:1/18:0), and Cer(d18:1/24:1) and lower concentrations of Cer(d18:1/24:0) were associated with a higher risk of fatal myocardial infarction [30].Our study found that Cer(d18:0/24:0) and Cer(d18:0/22:0) were significantly lower in T2DM patients with DR compared to those without DR, which suggests that different numbers of carbons and double bonds in ceramides might play differential roles in DR and CVD.The distinct ceramides and ceramide metabolites involved in metabolic regulation play unanticipated roles [31].Watt et al. discovered that circulating ceramides present in LDL particles were sufficient to induce insulin resistance in vitro and in vivo [32].However, how these two identified ceramides influence lipid metabolism in T2DM remains unclear and needs further exploration.Thus, disturbed Cer metabolism may contribute to dysfunction in DR, and therapeutic strategies to restore normal Cer metabolism might be an effective approach for treatment of DR.
In the discovery cohort, LPC(18:2) and LPC(16:0) were significantly higher in T2DM patients with DR.However, these two lipids were significantly lower in DR in the validation cohort.The previous findings point to a change in sphingolipid composition between control and T2DM [33].LPC is an inflammatory phospholipid and an important atherogenic substance in LDL that contributes to diabetic complications [34].Lipoprotein-associated phospholipase A2 (Lp-PLA2) plays a crucial role in diabetes-related retinal vasopermeability, a response mediated by LPC, and inhibiting Lp-PLA2 reduces diabetes-induced retinal vasopermeability [35].LPC O-16:0, LPC P-16:0, LPC O-18:0, and LPC 18:1 were all found to be inversely related to incident T2DM [36].The differences between the discovery and validation cohorts may be related to the populations studied, medications used, and stages of diabetic retinopathy [37].
There are some limitations of this study.First, only a Chinese ethnic group was selected, and future validation of our findings in other races or ethnic groups is warranted.Second, instead of chronic risk factors associated with the development of DR, some of the identified lipid markers might only represent temporary metabolic perturbations in this cross-sectional study.Third, the exact mechanism of DR development in patients with T2DM through which ceramide functions has not been explained.Therefore, more extensive preclinical and clinical studies are needed to clarify the mechanisms behind the potential effects of specific lipids.
Overall, the deregulation of sphingolipid metabolism in the diabetic retina appears to be a significant and seldom-studied element of DR pathophysiology.The precise mechanism underlying this disease is still unknown and requires further investigation.We showed the potential value of lipidomics research in understanding the pathophysiology of DR, and the results suggest that lipidomics profiling may be capable of identifying early-stage DR diagnostic indicators in high-risk Chinese populations.In addition, the findings from this study may help in the elucidation of new therapeutic targets for DR prevention and treatment.

Fig. 1
Fig. 1 Lipidome-derived markers identified from the discovery cohort.Lipidomic analysis identified fifteen candidate lipids of which serum levels were different between 42 T2DM patients with DR (DR group) and 42 T2DM patients without DR (NDR group) from the discovery cohort.(A) Mean peak intensity of lipids was analyzed after Log2 transformation of the data.(B) Fold change in DR/NDR was analyzed after Log2 transformation of the data.(C) Heatmap showing the distribution of lipid markers.Each row in the figure represents a different lipid, and each column represents a sample.Different colors indicate different intensities, and Log2 conversion was used for the data

Fig. 2
Fig. 2 The results of targeted lipidomics analysis in the validation cohort.For the validation cohort, the cases were 95 T2DM patients with DR (DR group), and the control subjects were 95 T2DM patients who had no DR (NDR group).(A) Peak area of lipids was analyzed after Log2 transformation of the data.(B) Fold change in DR/NDR was analyzed after Log2 transformation of the data.(C) The log2 conversion was used for the intensities of Cer(d18:0/24:0) and Cer(d18:0/22:0).All data are presented as the mean ± standard deviation (SD).Each symbol represents an individual participant.*p < 0.05, **p < 0.01, pairwise comparisons of change scores between the groups were evaluated by t test